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ABSTRACT 

We perform a systematic analysis of the effects of photometric redshift uncertainties on weak lensing 
tomography. We describe the photo- z distribution with a bias and Gaussian scatter that are allowed 
to vary arbitrarily between intervals of Sz = 0.1 in redshift. While the mere presence of bias and 
scatter does not substantially degrade dark energy information, uncertainties in both parameters do. 
For a fiducial next-generation survey each would need to be known to better than about 0.003 — 0.01 
in redshift for each interval in order to lead to less than a factor of 1.5 increase in the dark energy 
parameter errors. The more stringent requirement corresponds to a larger dark energy parameter 
space, when redshift variation in the equation of state of dark energy is allowed. Of order 10 4 — 10 5 
galaxies with spectroscopic redshifts fairly sampled from the source galaxy distribution will be needed 
to achieve this level of calibration. If the sample is composed of multiple galaxy types, a fair sample 
would be required for each. These requirements increase in stringency for more ambitious surveys; we 
quantify such scalings with a convenient fitting formula. No single aspect of a photometrically binned 
selection of galaxies such as their mean or median suffices, indicating that dark energy parameter 
determinations are sensitive to the shape and nature of outliers in the photo- z redshift distribution. 
Subject headings: cosmology - gravitational lensing, large-scale structure of the universe 



1. INTRODUCTION 

Weak gravitational lensing of galaxies by large 
scale structure is rapidly be coming one of the most 
powe rful cosmolo gical probes ijBartelmann fc Schneider! 
120011 iRefregierl 1200^ . Following the first detec- 
tions a few years ago iVan Wa erbeke et al.l 
Kaiser. Wilson fc Luppinol 120001 JBamn et alJ 



st uetec- 
aD J2000; 
lT l2000ai 



Kaiser, Wilson & luddhio zuuu: mco net au [zuuua 
Wittman et all 120001 IRefregier et al.l 12003^ . weak lens- 



ing has produced increasingly better constraints on the 



of mass fluctuations as (Hockstra. Yee &; Gladdcrs 2002 


Pen. van Waerbeke & Melliel .2002 


I iBrown et alJ 2003 


Jarvis et alJl2003l iPen et al.ll2003t 


Hevmans et alJl2005 


Van Waerbeke. Mcllier & Hockstra 


2005). While weak 



lensing is most sensitive to the amount and distribution 
of dark matter, it also has the potential to probe the 
dark energy throug h its effect on the growth of struc- 
tur e and distances lIHu fc Tegmarklll999t lHutererll2002t 
H^l2ll03llTakada fc ,Iainll2003HSong fc Knoxll2004t llshaH 
2005). Indeed, when combined with other cosmological 
probes, weak lensing data alre ady produce inter esting 
constraints on the dark energy l| Jarvis et al.ll2005^l . 

By utilizing source galaxy redshifts to study the growth 
of structure and the distance-redshift relation tomo- 
graphically, substantially more dark energy information 
can be recovered ()Hul ll999') . In fact future weak lensing 
surveys such as PanSTARRS 1 , Supe rnova/ Acceleration 
Probe rSNAP 2 : lATdering et alJl2004ft and Large Synop- 
tic Survey Telescope (LSST 3 ) are expected to impose 
constraints on dark energy that a re comparable to thos e 
from type la supernovae (see e.g. IRefregier et al.l l2004). 
In the more near term, the Canada-France-Hawaii Tele- 
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1 http:/ /pan-starrs. ifa.hawaii.edu 

2 http://snap.lbl.gov 

3 http://www.lsst.org 



scope Legacy Survey (CFHTLS 4 ) and the Dark Energy 
Survey 5 are expected to help bridge the gap between the 
current and ambitious future surveys. 

Powerful future surveys will require a much more 
stringent control of the systematics. Recent work 
has addressed systematic errors from t he computation 
of t he non-linear p o wer sp ectrum JVale fc White 
20031 IWhite fc Valel 120041: iHeitmann et all 1200 
Huterer fc Takadal 120051 lHagan. Ma fc Kravtsovl 12005ft . 



I 



baryonic cooling and press ure forces on t he distribution 
of la rge-scale structures ijWhitd 120041 iZhan fc Knoxl 
2004). ap proximations in inferrin g the shear from 
the maps UDodelson fc Zhand Iffl EHjj EM) ■ the 
presence of dust fVale et al.ll2004|) . Such studies have 
stimul ated work on how to improve the PSF reconstruc- 
tion ^aryjs_fc_Jain| E0Q4j), estimate shear from noisy 
maps (iBernst ein fc Jarvis! 2002; Hirata & Scliak 2003; 
iHoekstral l2004|) . and protec t against the small-sca le 
biases in the power spectrum l)Huterer fc Whitel f20051. 

In this work we consider the effect of errors in pho- 
tometric redshifts of source galaxies on weak lensing to- 
mography. Of course, the total number of galaxies, which 
is currently in the millions and might be in the billions 
with future surveys, is too large for spectroscopic mea- 
surements to be feasible. One therefore needs to rely on 
the photometric redshifts whose accuracy with the cur- 
rent state of the art in photometry, algorithms, galaxy 
classification etc, while presently adequate, may not be 
sufficient for future surveys which are expected to have 
very small statistical errors. Uncertain photometric red- 
shifts blur the tomographic bin divisions of source galax- 
ies. In the extreme case when photometric redshift errors 
are comparable to the width of the distribution itself, one 
completely loses tomographic information degrading the 

4 http:/ /www. cfht.hawaii.edu/Science/CFHLS 

5 http:/ /cosmology.astro. uiuc.edu/DES 



cosmological parameter accuracies by up to an order of 
magnitude. 

In this paper we study how the photometric redshift 
uncertainties affect cosmological parameter determina- 
tions. We construct an explicit mapping between the 
photometric and true redshifts, and parametrize it to al- 
low an arbitrary evolution of the bias and scatter between 
discrete redshift intervals. We then study how accurately 
the photometric redshifts need to be known a priori and, 
in particular, which details of the photometric redshift er- 
ror distribution are the main source of degeneracy with 
cosmological parameters. We hope that this study will 
help stimulate work on assessing and improving exist- 
ing algorithms for photometric redshift estimation (e.g 
ICunha et alll2005|) . 

The outline of the paper is as follows. In SJ21 we intro- 
duce the formalism and parametrizations of both cosmol- 
ogy and photometric redshift errors. We explore the loss 
of lensing information on the dark energy to photomet- 
ric redshift uncertainties in We show how this lost 
information is regained as we impose prior knowledge of 
the photometric redshift parameters in [f4] We discuss 
our results and conclude in <J21 

2. METHODOLOGY 

In this section, we discuss the modeling of the photo- 
metric redshift distribution. We then illustrate the flex- 
ibility of this description through two different fiducial 
models for the distribution. Finally we discuss its rela- 
tionship to lensing observables and the Fisher formalism 
for addressing its impact on parameter estimation. 

2.1. Photo-z Distribution 

Having only the photometric redshift ("photo-z") of 
the source galaxies at hand, the observer will necessarily 
bin the galaxies by their photometric redshifts z p h rather 
than true (spectroscopic) redshifts z. With a probability 
distribution p(z p \ l \z) in z p h at a given z, the true redshift 
distributions of the bins necessarily overlap. 

In general this distribution can vary arbitrarily with z. 
The true distribution of galaxies rii(z) that fall in the ith 



photo-z bin with 



"ph 



< Zph < 



^ph" 1 ^ becomes 



r z P h 

rii(z)= / dz ph n(z)p(z ph \z) . (1) 

n(z) = d 2 N/dzdfl is the overall galaxy redshift distribu- 
tion, and is chosen to have the form 



n(z) cc z a exp [— (z/z ) 



01 



(2) 



Unless otherwise stated we will adopt a = 2, (3 = 2 and 
fix zq such that median redshift is z mc d = 1. The total 
number of galaxies per steradian 



dzn(z) , 



(3) 



fixes the normalization, and we analogously define 

/>oo 

nf = dzrii(z) (4) 
Jo 

for the bins. 



By construction, the sum of the individual distribu- 
tions equals the total X)i n i( z ) = n ( z )- Therefore, re- 
gardless of how complicated the photo-z probability dis- 
tribution gets and hence the redshift distributions of the 
tomographic bins, the total distribution of galaxies n(z) 
is unchanged. 

This construction cleanly separates uncertainties due 
to the photometric redshifts of the individual survey 
galaxies characterized by p(z p h\z) from uncertainties in 
the redshift distribution of the underlying total distri- 
bution of galaxies n(z). We mainly consider the former 
in this work but comment on the latter in ^1 (see also 
IHuterer et a,l.ll2nM ITshak fc Hirata l l2nf)l. The ratio- 
nale is that even without any knowledge of the photo-z 's 
of the survey galaxies themselves, one can at least bin 
all of the galaxies together assuming that the underlying 
redshift distribution or selection function of the survey 
is known. In practice, this means that one must obtain 
information about the underlying distribution from an 
independent source (say, another survey through a study 
of the luminosity function) or from a fair subsample of 
survey galaxies with spectroscopic redshifts. 
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Fig. 1. — Parametrization of the photo-z distribution and two 
illustrative fiducial models. The distribution spreads galaxies at 
a given redshift z into a distribution in z p h characterized here by 
a bias zt,i as and a scatter a z whose evolution is parametrized by 
interpolating their values at N pz r edshifts z M . In both models 
2 bias = 0, whereas a z , given in i|2.2l is illustrated here for model I 
(shaded region) and model II (dashed lines) as lcr bands. Galaxies 
binned according to their photometric redshifts (-/V<jj v horizontal 
bands) then have overlapping redshift distributions determined by 
the 27V pz photo-z parameters. 



2.2. Photo-z Models 

Any photo-z model may be described by providing 
a function for the distribution of photometric redshifts 
given the true redshift, p(z p h|z). For the purposes of 
this paper we take the simplifying assumption that this 
function is a Gaussian at each redshift, i.e. 



P(z p h\z) 



1 



2na 



■ exp 



(-2 ^ph ^bias) 



2g 7 , 



(5) 



2 



However, we allow the bias Zbias(z) and scatter u z (z) to 
be arbitrary functions of redshift. The rcdshift distribu- 
tion of the tomographic bins defined by equation Q can 
then be written as 



rn{z) = \n(z) [erf(x i+ i) - erf(a; i )], 



with 



_/ (») 
x l = (z ph 



Z + Zbia 



2ct 2 



(6) 



(7) 



where erf (x) is the error function. 

The Gaussian assumption is not as restrictive as it 
might naively seem. By allowing the bias and scatter 
to be arbitrary functions of redshift one can obtain arbi- 
trarily complex redshift distributions in the tomographic 
bins through equation (Q. In fact, the mapping is in 
principle completely general for finite bins and a smooth 
underlying distribution. Galaxies in a finite range of red- 
shift over which the distribution is nearly constant can 
then be mapped to any z p h- 

In practice we will represent the free functions Zbi as (z) 
and a z (z) with a discrete set of N pz photo-z parameters. 
They represent the values of the functions at z M which are 
equally spaced from z — to 3. To evaluate the functions 
at an arbitrary redshift, we take a linear interpolation of 
the discrete parameters in redshift. 

While a finite N pz does restrict the form of the dis- 
tribution, it still allows radically different redshift distri- 
butions given the same tomographic bins. For example, 
consider two different photo-z models 

• Model I: z hias ( z ) = 0; °z{z) = 0.05(1 + z). 

• Model II: z hias (z) = 0; a z (z) = 0.2 for z < 1.0 and 
a z (z) = 0.5 for z > 1.0. 

The distribution p(z p h|z) is illustrated in Figure ^ f° r 
N pz = 31 through the la scatter region. The resulting 
redshift distributions for N^w = 5 tomographic bins are 
shown in Figure^ These specific choices of iVp Z and N^w 
are motivated in 

Model II demonstrates that sharp changes in the Gaus- 
sian photometric parameters can map neighboring galax- 
ies in redshift to quite different tomographic bins. The 
redshift distributions of the bins can thus have features 
that are sharper than the assumed scatter. Addition- 
ally, photo-z degeneracies that take two distinct spectro- 
scopic redshift ranges into a single photometric redshift 
and lead to bimodality in the binned distribution can 
be modeled by a large Zbias- Finally, galaxy types with 
different photo-z distribution at a given redshift can be 
approximated by discontinuous jumps between infinites- 
imally spaced redshift bins. Such considerations, while 
potentially important are beyond the scope of this work. 

In summary allowing our set of parameters to freely 
vary one can access a wide range of tomographic redshift 
distributions. Uncertainty in these parameters will then 
cause uncertainties in tomographic dark energy determi- 
nations. 

2.3. Lensing Observables 

The convergence power spectrum at a fixed multipole 
I and for the «th and j th tomographic bin (£) is given 
hv llKaiseHI1992irm9l 

r°° h(z) 

nfnfP«(e)=J^ dzWiWWjWp^Pfaz), (8) 
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Fig. 2. — Source galaxy redshift distribution n(z). Top panel: 
photo- z model I. Lower panel: photo-z model II. The solid curve is 
the overall galaxy distribution defined in equation The other 
curves are the true (spectroscopic) distributions that correspond 
to the sharp divisions in photo-z space (denoted by dotted vertical 
lines). 



where H(z) is the Hubble parameter, and D(z) is 
the angular diameter distance in comoving coordinates. 
P(k(;, z) is the three-dimensional matter power spectrum 
and kg = £/D(z) is the wavenumber that projects onto 
the multipole t at redshift z. The weights W are given 

by 



, , 3 H2D(z) , 

Wi(z) = -n m ^^(i + z) 



x / dz'ni(z') 



Dls{z,z') 
D(z') 



(9) 



where D L s(z,z') is the angular diameter distance be- 
tween the two redshifts. The power s pectrum is com- 
puted from the transfer function of lEisenstein fc Hul 
(1999) with dark energy modification s fromlHul (12 002). 
and th e non-linear fitting function of iPeacock fc Doddsl 
(.1993). 

With tomographic binning, the number weighted 
power spectrum nfn^P^ and not is the fundamental 
observable. Even given photometric redshift uncertain- 
ties in the binning it is always possible to recover the 
total weighted power spectrum l)Hul ll999) 



(n A ) 2 P K = J2 nfnfP* 



(10) 



since the weighting is based on the observed nf~. By 
treating nfn^Pj* as the observable one guarantees that 
the addition of photo-z estimates for the individual 
galaxies can only add information. This would not be 
true if P£j were taken as the only observable quan- 
tity. Given that changes in photo-z parameters induce 
changes in nf , the binned power spectra PX do not con- 
tain enough information to weight the power spectra and 
recover the total P K . 

That the binned angular number densities are ob- 
served quantities also implies that there is additional di- 
rect information on the photo-z parameters that does 
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not depend on shear measurements. For example, a 
high fraction of galaxies in bins with z p h larger than the 
median redshift would imply a large photo-z bias. We 
choose not to consider this sort of information since it is 
not directly related to lensing. Furthermore for the small 
changes in nf that we will typically be considering, the 
sample variance between the observed and that pre- 
dicted by the underlying redshift distr ibution and the 
photo -z parameters cannot be ignored l|Hu fc Kravtsovl 
12001 . Therefore we will consider the number weighted 
power spectra nfn^P^ as the fundamental lensing ob- 
servables. 

2.4. Fisher Matrix 

The Fisher matrix quantifies the information contained 
in the lensing observables 

O a = i{i ^ )/2+j {l)=ntnfP? j {l), (i>j) (11) 

on a set of cosmological and photo-z parameters p M . Un- 
der the approximation that the shear fields are Gaussian 
out to £ max , the Fisher matrix is given by 



F. v = g X (2^+l)/ sky ^^[C- 1 ] afc a ° b 



1=2 



(12) 



so that the errors on the parameters are given by Ap^ = 

Given shot and Gaussian sample variance, the covari- 
ance matrix of the observables becomes 

Cab = nfnfnint (F^Pj? + P^Pjk) , (13) 

where a = i(i- l)/2 + j, b = k(k - l)/2 + I. The total 
power spectrum is given by 



ptot 

ij 



I hi 



(14) 



where 7; n t is the rms shear error per galaxy per compo- 
nent contributed by intrinsic ellipticity and measurement 
error. For illustrative purposes we will use f max = 3000, 
/sky corresponding to 4000 sq. deg, n A corresponding to 
55 galaxies/arcmin 2 and 7i n t = 0.4. The value of £ max is 
motivated by simulations which find substantial devia- 
tions from Gaussianit y and weak lensing approximations 
on ar cminute scales l)White fc Hul 120001 iVale k, White I 
2003). For our choice of noise parameters above, the re- 
sults are insensitive to ^ max since the measurements are 
noise dominated on those scales. 

For the cosmological parameters, we consider four pa- 
rameters that affect the matter power spectrum: the 
physical matter density il m h 2 (— 0.14), physical baryon 
density f2f,/i 2 (= 0.024), tilt n s (= 1), and the amplitud e 
5 C (= 5.07 x 10~ 5 ; or A = 0.933 ISoergel et alJ lt2003h l. 
Values in parentheses are those of the fiducial model. 
Unless otherwise stated, we shall take priors on these 
four parameters of AlnO m /i 2 = A In ilf, h 2 = A In 8^ = 
An s = 0.05. These priors represent only a modest im- 
provement over current determinations. Our results on 
the relative degradation in constraints caused by photo- 
z errors are insensitive to reasonable variations in this 
choice. 

To these four cosmological parameters, we add either 
two or three dark energy parameters: the dark energy 



density f2oE(= 0.73), its equation of state today wq = 
Pde/pde|z=o(= —1) and optionally its derivative w a = 
—dw/da\ z= Q(~ 0) assuming a linear evolution with the 
scale factor w = wq + (1 — a)w a . 

Note that throughout this paper, our notational con- 
vention is latin indices for tomographic bins and greek 
indices for parameters. 
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Fig. 3. — Relative errors in dark energy parameters as a func- 
tion of the number of tomographic divisions iV^jiv Solid lines cor- 
respond to w a for the set {aiciOai^DE}; dashed lines to wq for 
{wo, Hde}- Both lines assume that all photo-2 parameters are per- 
fectly known (i.e. fixed). Note that the results converge at smaller 
JVjiv for a smaller dark energy space, and that N^w = 5 is more 
than sufficient in either case. The points correspond to the same 
cases, but now with 27V pz = 62 photo-z parameters marginalized 
(Photo-z priors of Aa z = Azi, iaa = 10 are applied). Here essen- 
tially all tomographic information is lost so that the errors are 
comparable to those of N^i v = 1 or no tomography. 



3. DARK ENERGY INFORMATION LOSS 

In this section we consider the nature of the tomo- 
graphic information on the dark energy and its loss to 
photo-z uncertainties. We establish the maximal infor- 
mation that can be gained through tomographic redshift 
divisions for a given dark energy parametrization. We 
then determine the number of photo-z degrees of freedom 
that would be required to lose this information. This loss 
of information is caused by a degeneracy between cosmo- 
logical and photo-z parameters. We explicitly construct 
an example of this degeneracy as both an illustration and 
test of our statistical methodology. 

3.1. Maximal Information and N,a v 

For any given choice of dark energy parameterization, 
the information contained in lensing will sat urate wit h 
some finite number of tomographic bins A^jiv <|Hul [l999). 
Since the broad lensing kernel of equation 10 makes the 
shear for neighboring source redshifts highly correlated, 
most of the information is contained in a few coarse bins. 
The exact number depends on the type of information 
that is to be extracted. Roughly speaking, the number 
of bins should exceed the number of dark energy param- 
eters. 

Figure|3 (lines) quantifies this expectation through im- 
provement in the errors on dark energy parameters as a 
function of iVdiv for Model I. For a 2 parameter dark en- 
ergy space {wo,^de}, ^Vdiv = 3 divisions equally spaced 
from z = to z = 3 are enough for the improvements 
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TABLE 1 

Baseline constraints on dark energy parameters 



Photo-2 Model 


Parameters 


cr(fi DE ) 


a(w ) 




I 


{n DE ,ui } 


0.0062 


0.061 






{0 D Ei Wo, 1tl a } 


0.024 


0.25 


0.69 


11 


{n DE ,i«o} 


0.0073 


0.070 








0.034 


0.36 


0.96 



in wo to saturate. For a 3 parameter dark energy space 
{wq, w a , ^de}, Ndiv = 4 divisions are sufficient for w a . 
Note that iVdj v — 1 corresponds to no tomography or 
no photo-z information on the individual galaxies. The 
dark energy parameters that are not shown in Figure 
behave similarly. In what follows we conservatively adopt 
Ndiv = 5 as sufficient to extract the dark energy infor- 
mation. With Ndiv = 5 and photo-z parameters fixed, 
the constraints on dark energy parameters are shown in 
Tabled 

Note that improvements relative to the no-tomography 
case are more significant in the larger parameter space. 
This is due to the fact that wq is nearly degenerate with 
w a since lensing mainly constrains w(z) at some inter- 
mediate redshift (see below). Even the small amount 
of information in the fine-binned tomography assists the 
breaking of the degeneracy. 

3.2. Maximal Degradation and N pz 

Next we choose the number of photo-z parameters N pz 
that describe each of the functions Zbi as (z) an d & z (z). We 
seek to allow enough freedom in the photo-z parameters 
so that in the absence of prior information on their val- 
ues all of the tomographic information is lost. Because 
the limit of no tomographic information corresponds to 
Ndiv = 1, we have a quantitative means of assessing the 
minimal N pz . When N pz becomes large enough, the vari- 
ations in redshift, which act on the characteristic scale 
of Sz = 3.0/ (N pz — 1), are rapid enough that they do not 
mimic any variation in cosmological parameters. 

Figure 0] shows the degradation in the errors on w a for 
the cases of N pz = 6, 11, 21, 31 and 61 as a function 
of the prior on the photo-z parameters. Results for wq 
are similar. To compare priors for different N pz values, 
we have here rescaled the individual parameter priors by 
y/N pz /31 so as to be equal for a fixed redshift interval 
Sz = 0.1. The results have converged with N pz > 21. To 
be conservative, in the rest of this paper we use N pz = 31, 
or a total of 62 photo-z parameters. 

The impact of this choice of N pz = 31 as a function 
of TVdiv for dark energy parameters is shown in Figure |3 
(points). For all A*di V , these constraints match those with 
no tomographic binning very well, showing that without 
prior information on the photo-z parameters all tomo- 
graphic information has been effectively destroyed and 
we recover the case with a single redshift division. The 
small discrepancy comes from the inadequacies in the 
Fisher matrix one of which is the local approximation 
to the parameter errors as we shall discuss in the next 
section. 

3.3. Photo-z - Dark Energy Degeneracy 




Fig. 4. — Error degradations in w a (that is, errors in w a relative 
to the error with perfect knowledge of photo-z parameters) as a 
function of the photo-2 prior. The photo-z priors are rescaled by a 
factor of \J iVp Z /31 so that they reflect constraints per Sz = 0.1 in- 
dependently of iVp Z . Different lines from top to bottom correspond 
to different iV pz : 61 (short dashed line), 31 (solid line), 21 (dash 
3-dotted line), 11 (dash dotted line) and 6 (long dashed line). Note 



that the results have converged with JV pz 
just to be conservative. 



> 21: we use N u 



31 



With a sufficient number of unknown photo-z parame- 
ters 2N pz > 62, the Fisher matrix results of the previous 
section imply that dark energy information in tomogra- 
phy is completely lost. This fact implies that there is a 
nearly perfect degeneracy between photo-z, dark energy 
and other cosmological parameters. Here we examine 
that aspect of the degeneracy that involves the photo-z 
and dark energy parameters only. This degeneracy alone 
suffices to destroy most of the tomographic information 
and will remain even if the other cosmological parameters 
are perfectly measured from other sources. 

Constructed from parameter derivatives, the Fisher 
matrix is a local expression of the degeneracy in pa- 
rameter space. Because the Fisher matrix results im- 
ply that the degeneracy persists to large changes in the 
dark energy parameters, it is important to assess the ex- 
tent of the degeneracy more directly and test the validity 
of the Fisher approximation. If the degeneracy relation 
"curves" in parameter space, the Fisher approximation 
will only find the local tangent. 

We start by identifying this local tangent with the 
Fisher matrix. To isolate the degeneracy between dark 
energy and photo-z parameters, we eliminate the other 
cosmological parameters, formally by adding strong pri- 
ors to the Fisher matrix. For numerical reasons we 
also add a weak prior on photo-z parameters (Azbi as = 
A<7 2 = 1 ) to control numerical errors from the nearly sin- 
gular Fisher matrix. Of the eigenvectors of this Fisher 
matrix that involve the dark energy, those with the small- 
est eigenvalues will be responsible for most of the photo- 
z dark energy degeneracy. We find that a single linear 
combination of parameters (dark energy and photo-z) 
contributes most (~ 98%) of the errors in dark energy 
parameters wq and w a . Thus the degeneracy is essen- 
tially one dimensional in the multi-dimensional parame- 
ter space. Let us call this direction - or the eigenvector 
of the Fisher matrix - e w . 

The true extent of the degeneracy is quantified by the 
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change in \ 2 between the fiducial model p^ and a trial 
model p M 



£=2 



ab 



x [CT 1 ].,, [Ot^Pp) - O b {£;p^)} . (15) 

If the Fisher matrix approximation were valid out to say 
la along the degeneracy then the trial model p^ = p^ + 
a w e w , where a~ 2 is the eigenvalue corresponding to e w , 
would be separated by 



CF—^Fe w = l (16) 

due to the orthonormality of the eigenvectors. In prac- 
tice Ax 2 rnc — 857 for this extrapolation indicating a cur- 
vature in the degeneracy direction. In other words, p^ 
and p^ are highly distinguishable models in spite of the 
Fisher prediction that they are indistinguishable. 

Even given curvature in the degeneracy direction, the 
Fisher approximation remains useful if it accurately pre- 
dicts the extent of the degeneracy. This is especially true 
if the curvature lies mainly in the photo-z nuisance pa- 
rameters which exist only to be marginalized. To assess 
the extent of the degeneracy, we use the Fisher matrix as 
a local approximation of the degeneracy with the follow- 
ing procedure. Starting at the fiducial model, calculate 
the Fisher matrix and find e w as defined above, then take 
a small step along e w direction. Now calculate the Fisher 
matrix at the new point, find the new e M , take another 
small step along this new direction. Repeat. The small- 
ness of the steps is controlled such that the change of \ 2 
between steps agrees within 10% with that predicted by 
Fisher matrix. This 10% error may add up to a much 
bigger percentage after a few steps. To make sure that 
we stay on the degeneracy direction, local minimum of 
Axtruc is searched after a few steps out. 




Fig. 5. — Equations of state of dark energy w(z) that are de- 
generate with photo-z parameters. Solid lines show a series of 
degenerate models stepping out by Auio = 0.2 and ending at a 
model with Aiuo = 1.14 which deviates from the fiducial model at 
+lcr for the fiducial survey. The tight correlation between wq and 
w a along the degeneracy direction results in the tight "waist" or 
pivot where w(z) remains well determined. The dashed line shows 
that the +lcr degenerate model as predicted by the Fisher matrix 
is in good agreement with the true degeneracy even out to large 
Aw . 

From this construction we find that the extent of the 
degeneracy in w is accurately predicted by the Fisher 



matrix. Figure [3] shows that the model with A^p = 1 
(thick dashed line) is almost identical to the model with 
AXtruc = 1 (thick solid line) in w{z). In Figure we 
also show intermediate models along the degeneracy with 
Axtruc < !• That they all pass through essentially a sin- 
gle point in w(z) space is another indication that the de- 
generate direction lies almost entirely along a specific lin- 
ear combination of dark energy parameters as predicted 
by the Fisher matrix. The curvature in parameter space 
mainly involves the photo-z parameters. 

The redshift at which these curves intersect is z « 0.7 
and at this redshift measurements of w are essentially 
immune to photometric redshift errors. This immunity 
reflects the fact that even without tomography lensing 
can constrain the equation of state at some effective red- 
shift. With two parameters to describe w(z), there is 
only one remaining linear combination to be affected by 
photometric redshifts. With a more general parameteri- 
zation of w(z) we expect that there will be multiple de- 
generate directions with roughly the same single aspect 
of w(z) preserved. 

Furthermore, the two estimates (true and Fisher) agree 
on the amplitude of the photo-z parameter variation 
along the degenerate direction. For example, as shown 
in Figure|Sl at a point along the e w direction which is la 
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Fig. 6. — Comparison of the fractional differences of the redshift 

distribution Srij/n in each of the i = 1, , 5 tomographic bins that 

establishes the la degeneracy with dark energy parameters of an 
extent Auio ~ 1- Also shown are the changes in the two sets of 
discrete photo- 2 parameters <5pz = (<5<t z , <5zbias) that cause these 
differences, connected by lines to guide the eye. Top panel (F): 
Fisher approximation for the lcr degeneracy. Bottom panel (T): 
true degeneracy. 



away from the fiducial model, the Fisher matrix indicates 
that the photo-z parameters changed by <5zbi as < 0.06 
and Sa z < 0.07, while the actual bounds on the vari- 
ations are Sz^ is , s < 0.04 and Sa z < 0.07. Note that 
these changes are fairly small and imply that subtle vari- 
ations in the redshift distributions for the tomographic 
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bins are responsible for a degeneracy that degrades errors 
in wq and w a by an order of magnitude. Figure El also 
shows that the difference between these distributions for 
Axlrue = 1- We expect that with a change in the photo-z 
model, the specific photo-z variations that establish this 
degeneracy will change but that a strong degeneracy will 
remain. 

In summary, we find that the Fisher matrix is an ade- 
quate tool for assessing the existence and extent of degen- 
eracies between photo-z and dark energy parameters. It 
should not however be used to infer the specific changes 
in the photo-z parameters that establish the degeneracy. 

4. PHOTO-Z INFORMATION RECOVERY 

In the previous section, we established the existence 
of a degeneracy between photo-z parameters and dark 
energy parameters and tested the validity of the local 
Fisher matrix approximation to this degeneracy. In this 
section, we use the Fisher matrix formalism to investi- 
gate the extent to which prior information on the photo-z 
distributions help recover the tomographic dark energy 
information. We assume 2N pz = 62 photo-z parameters 



and iVdi, 
(see Q. 



5 tomographic bins through out this section 




Fig. 7. — Error degradations for constant wo (left panel) and 
those for w a when both wq and w a are varied (right panel), as a 
function of photo-z parameter priors. Here the degradations are 
defined as actual errors relative to errors that assume the photo-z 
parameters to be perfectly known. Priors on the photo-z param- 
eters Zbias(Z(j) a n d 0"z(z/j) are shown on the x-axis and y-axis re- 
spectively. Fiducial photo-z model I is assumed and the photo-z 
parameter spacing in redshift is <5z = 0.1. 



4.1. Photo-z Priors 

We now explore the effect on dark energy parameter 
constraints of priors on each of the photo-z parameters 
Zbias(zp) and a z (z fl ). For simplicity, we begin by ap- 
plying a redshift independent prior on the parameters. 
In practice parameters controlling the distributions well 
above and well below the median redshift require weaker 
priors. We will discuss this point in i j4.3l 

In Figure the left panel shows the error degrada- 
tion in wq assuming the {wq, ^de} parametrization (left 
panel) and w a assuming the {wq, w a , ^de} parametriza- 
tion (right panel). For reference the baseline errors for 
the fiducial survey are listed in Tabled 



As in the previous section where other cosmologi- 
cal parameters were artificially fixed, we find that the 
larger dark energy parameter space is more susceptible 
to photo-z errors. For example, for the extreme case of 
no photo-z information, i.e. very weak priors on both 
bias and scatter parameters, dark energy parameters of 
the {wo^de} parameterization are degraded by about 
a factor of two while those of {wq, w a , ^de} parameteri- 
zation are degraded by about a factor of ten. 

In the more relevant case where we demand that the 
dark energy error degradation be no larger than 1.5, 
the requirement per photo-z parameter is about 0.01 for 
{wq, Ode} and 0.003 for the {wg, w a , ^de} parametriza- 
tion. Figure [7| also shows that both bias and scatter pa- 
rameters are important, and that knowledge of the bias 
is only slightly more important than that of the scatter. 
Furthermore, dark energy parameters that are not shown 
in Figure have very similar requirements to those plot- 
ted in either parametrization. 

4.2. Dependence on Fiducial Model 

The results above are for a specific choice of the fidu- 
cial model for the photo-z distribution and survey. To 
explore the dependence on the former we take the very 
different photo-z Model II, where the scatter is substan- 
tially larger and jumps discontinuously in redshift. Even 
so the requirements on the photo-z parameters are very 
similar; see Figure [5] In particular, within the interest- 
ing regime where the photo-z prior is smaller than unity 
and the degradation in dark energy parameter errors is a 
factor of a few or lower, the two models agree very well. 

On the other hand, photo-z requirements do depend 
on the parameters / s ky, 7int and n A that determine the 
level of sample and noise variance in the survey. The 
trend is that the more ambitious the survey, the more 
stringent the requirements on photo-z parameters. The 
scaling for the priors on photo-z parameters can roughly 
be described as 

A P(/sky,7int,n A ,d) 

Ap(0.1,0.4,55,1.5) " 




gives the photo-z pri- 
Here g(d) 



where Ap = Az h - ms = 
ors and n A is in units of arcmin 
Ap(0.1, 0.4, 55, d) scales the prior requirement to alter- 
nate levels of degradation d; it is shown in Figure as 
d(g) and is only weakly dependent on the fiducial photo- 
z model. With the best fit C = 0.6 this scaling is good 
up to a factor of 2 for any reasonable set of survey pa- 
rameters. 

Finally the photo-z precision requirement is not very 
sensitive to z me d, the median source redshift of the sur- 
vey. For 0.68 < z moc i < 1.3, Ap varies by less than 40%. 

4.3. Training Set Size 

The priors on the photo-z parameters ultimately re- 
quire a training set of galaxies with measured spectro- 
scopic redshifts. Operationally suppose that a photo-z 
training set has N£ spectroscopic redshifts per redshift 
interval determined by N pz (here Sz = 0.1). 

Given a Gaussian distribution for the photo-z distribu- 
tion and a fair sample of spectroscopic galaxies selected 
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Fig. 8. — Comparison of the photo-z requirements for our two 
fiducial photo-z models. The solid line shows the degradations 
for model I while the dashed line corresponds to model II. The 
fiducial errors in w a when the photo-z's are perfectly known are 
(70(11)0) = 0.69 (model I) and (jQ(w a ) = 0.96 (model II). 



from this distribution, the training set would indepen- 
dently determine the bias and scatter to 



AZbias(^) =Cr z (z fi )yl/Ntpec , 

Aa z {z^=a z {z^^2/Nt pcc . (18) 

For a fixed dark energy degradation, N^ pcc depends on 
two things: the fiducial a z and the required prior as 
scaled from equation (|17|l . Since the photo-z prior re- 
quirement is roughly independent of a z as shown in Fig- 
ure|Sl the larger the scatter, the larger the required train- 
ing set. Note that A/gpec (= J^p. ^s P cc ) i s robust to 
changes in the number of photo-z parameters or Sz. For 
example a binning of Sz = 0.05 would imply twice as 
many photo-z parameters that would need to be con- 
strained a factor of V2 less well, yielding the same re- 
quirement on N spcc . 

For determining the redshift extent of the training set, 
it is important to go beyond our simple redshift indepen- 
dent prior. FigurelTUI shows the cumulative N spcc (> z) 
required for 1.5 degradation in dark energy. Notice that 
this flat prior assumption would require a substantial 
number of galaxies across the whole redshift range (10 5 ), 
including 8 x 10 4 galaxies above z = 1.5. This number 
is artificially high since the actual requirements on the 
prior fall sharply away from the median redshift of the 
distribution. 

To illustrate the difference, we constructed a weighted 
template of how the bias and scatter priors vary with 
redshift to produce a fixed degradation in dark energy 
parameters. We choose a simple power law of n.; which is 
the number of galaxies in each redshift interval (Sz = 0.1 
for the fiducial choice of N pz = 31) corresponding to the 
photo-z parameters. At z < 1.2, the power is chosen as 
— 1. To account for the difficulty in measuring redshifts 
at z > 1.5 from optical bands, we steepen the index 
to —3 for z > 1.2. Figure^ compares the flat prior to 
the weighted one. The requirement for the weighted one 
drops to a total of 4 x 10 4 but more importantly only 300 
at z > 1.5. 

For dark energy degradations other than 1.5, we pro- 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Fig. 9. — Photo-z prior templates of how the bias and scatter pri- 
ors vary with redshift to produce a fixed error degradation in dark 
energy parameters. The degradation on w a is 1.5. The straight 
line (1; dashed) is the flat template used in FigurelTl The curved 
line (2; solid) is the weighted template constructed according to the 
galaxy number density. At z < 1.2, the curve follows n~ while at 
z > 1.2, the curve follows n~ 3 . Here fii is the number of galaxies 
in each redshift interval Sz = 0.1. 



vide in FigurelTTlthe ratio of A^pe,, for an arbitrary dark 
energy degradation to that of 1.5 degradation. In order 
to find out the N spcc requirement for any dark energy 
degradation, all one needs to do is to look up the ratio 
in FieurelTTl and multiply it by the -/V spC c in FigureTTTil 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



z 

Fig. 10. — Cumulative requirement of N BpC c for 1.5 dark energy 
error degradation. The corresponding photo-z prior templates are 
shown in FigurelTfl with (1; dashed) as the flat prior and (2; solid) 
as the weighted prior. 

We have tested that the scaling relation to surveys 
with different fiducial parameters of equation (|17f) works 
equally well for both the flat and weighted priors. Using 
the scaling relation, the requirement of N spcc could be 
scaled to different survey easily. 

4.4. Mean vs. Median 

We now estimate the amount of information that 
comes from knowing the median or mean of the redshift 
distribution of source galaxies in each tomographic bin. 
Note that this is distinct from priors on the photo-z bias 
or mean photo-z at a given redshift. This question is 
interesting in its own right, but also because other work, 
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12 3 4 

d = <j(w ) / a„(w ) 



Fig. 11. — The ratio of required N Bpcc for an arbitrary w a degra- 
dation relative to that of 1.5 degradation for the two prior tem- 
plates of Figurdi^l 
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Fig. 12. — The effect of the mean or median of tomography bins. 
Vertical axis is dark energy error degradation. The solid line and 
the dashed line are for the case of mean and median respectively. 
Photo- z priors of Aa z = Az^ias = 1 are assumed. 



parallel to this l|Huterer et alJ l2005j) has parametrized 
the photo-z uncertainty by the centroids of the tomo- 
graphic bins, which when varied shift the overall distri- 
bution of the corresponding tomographic bin. While one 
intuitively expects that the centroid of the photo-z bin 
(or, more generally its mean) carries the most informa- 
tion, we now have the tools to precisely examine the rel- 
ative contribution of the mean relative to that of the 
higher moments. 

Figure 1121 shows degradations of the error in w a as a 
function of priors on the mean or median of the tomo- 
graphic bins. All of the photo-z parameters are given a 
weak prior of unity for numerical stability. A prior of 
10~ 3 on the mean is enough to render the mean essen- 
tially precisely known. But the dark energy degradation 
is still a factor of 2 even with perfect mean measurements 
showing that there remains information lost to the higher 
moments of the distribution. 

Similarly, while the mean of the tomographic distribu- 
tion does carry the majority of the information, the me- 
dian carries significantly less (see Figure H^I . The reason 
is that the mean has extra information about the tails of 
the redshift space distribution while the median does not. 



This sensitivity to the tails will make obtaining precise 
measurements of the mean difficult. One still requires a 
fair sample from each of the tomographic redshift bins 
extending to high redshift. In the end, mean priors re- 
quire a similar number of training set galaxies as above. 

5. DISCUSSION 

We have performed a systematic study of the effects of 
imperfect photometric redshifts on weak lensing tomog- 
raphy. Describing the photo-z distribution with a bias 
and scatter that can vary arbitrarily between redshift in- 
tervals of 5z — 0.1, we studied the degeneracies between 
photo-z and dark energy parameters, as well as the re- 
sulting degradations in dark energy parameter errors. 

Not surprisingly, we find that there exist significant de- 
generacies between the dark energy and photo-z param- 
eters. Assuming that the overall distribution of galaxies 
n(z) is independently known and the photometric red- 
shifts are used only for tomographic subdivision, we find 
that larger dark energy spaces suffer more degeneracy 
with photo-z than the smaller ones. 

Without any information on photo-z parameters, one 
recovers the no tomography case where errors on fidu- 
cial parameters are a factor of two times worse (for the 
{woj^de} parametrization) or ten times worse (for the 
{wq, w a , ^de} parametrization) than those for the 10-bin 
tomography case with perfect photo-z's. 

For the fiducial survey, in order to have less than a fac- 
tor of 1.5 degradation in dark energy parameter errors, 
the uncertainties of photo-z parameters Zbi as arid a z (de- 
fined in the redshift interval Sz = 0.1) should each be 
controlled to better than 0.003-0.01, depending again on 
the size of dark energy parameter space. We provide a 
convenient approximation for scaling these requirements 
to different surveys. Importantly, no single number such 
as the mean or median of galaxies in the tomographic 
bin captures all of the effect of photo-z errors. That the 
mean captures more of the information than the median 
indicates that the dark energy information is sensitive to 
the tails of the distribution. 

In order to achieve less than a factor of 1.5 degradation 
in the evolution of the equation of state, a training set 
of a few times 10 4 galaxies with spectroscopic redshifts 
is required. Again, one can easily rescale the number of 
galaxy requirement to different surveys using our scaling 
relation. 

There are several caveats to our assessment that merit 
future study. Although our parametrization can handle 
photo-z degeneracies, for example from multiple galaxy 
types, that lead to bimodality and catastrophic errors, we 
have limited our study to fiducial models around which 
their effects are small. Such effects will increase the re- 
quired number of training set galaxies. Moreover, we 
have assumed that the parent redshift distribution of the 
survey is known and that photometric redshifts are only 
employed to subdivide the galaxy sample for tomogra- 
phy. Uncertainties in the parent distribution can further 
degrade dark energy determinations. 

On the other hand, uncertainties in the parent dis- 
tribution are also constrained by the training set. If 
one assumes that n(z) is a smooth function that is 
parametrized by relatively few parameters, uncertain- 
ties in the parent distribution should be smaller than 
those of the tomographic bins. For illustrative purposes, 
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if we parametrize n(z) with the three parameters of equa- 
tion J2J), we find that the constraint on n(z) from N^ pcct 
is good enough to have dark energy parameter errors dif- 
fer by less than 10% from the case where n(z) is fixed. 

Given the current state-of-the-art of photo-z algo- 
rithms as well as expected improvements with multi- 
wavelength observations of all source galaxies, prospects 
for sufficiently accurate determination of photometric 
rcdshifts are bright. Nevertheless, it will be an important 
and challenging task to achieve good control of the photo- 
z accuracy for the specific types of galaxies selected in 



lensing surveys, and then propagate the remaining photo- 
z errors into the final cosmological constraints. 
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